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A Monte Carlo simulation method has been developed for 
studying nuclear-electromagnetic cascades initiated by high 
energy nuclei with Z > 3 incident on heavy absorbers. The 
calculations are based on a cascade model which was first 
adjusted until it agreed with measurements made with protons 
at an accelerator. Modifications of the model used for pro- 
tons include the Incorporation of the probabilities for frag- 
mentation of heavy nuclei into lighter nuclei, alpha particles 
and nucleons. Mean values and fluctuations of the equivalent 
numbers of particles in the cascades at various depths in an 
iron absorber are presented for protons, carbon, and iron 
nuclei over the 30 to 300 GeV/nucleon energy range. 

K Introduction. To be presented here are some of the first 

results from a Monto Carlo study of the nuclear-electromagnetic cas- 
cade development initiated by high energy nuclei with Z > 3 inci- 
dent on heavy absorbers. The present results are for C and Fe 
nuclei incident on an iron absorber. The absorber, in keeping 
the design of ionization spectrometers, conslsbs of alternate 
of iron and plastic scintillator. The iron layers are each 2 ra- 
diation lengths (r.l.) thick, and the scintillators, used as the 
cascade sampling layers, are each 1/4 in. thick. The cross-sectl 
al area used for the absorber is 50 x 50 cm"^. 

The cascade model used in the calculations is an extension 
of a model originally devised for studying the cascade development 
initiated by protons in a small iron absorber ionization spectro- 
meter, which was exposed both to cosmic rays in a series of balloon 
flights (Schmidt, et al. I969) and to artificially accelerated pro- 
tons (Jones, et al. I969) from the Brookhaven Alternating Gradient 
Synchrotron (AGS). The model for protons was first adjusted until 
it agreed with the AGS measurements. Calculations were then carried 
out for higher incident energies and for absorbers other than iron 
(Jones, 1970). Later, the model was also extended to predict the 
cascade development initiated by alpha particles (Jones, 1971). 

The present work being done on the cascade development initiated 
by nuclei with Z > 3 has been undertaken because of the Interest 
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in using ionization spectrometers for studying the energy spectra 
of these nuclei in the primary cosmic rays. 

:?■ Calculations. Space does not permit a detailed description 

ol^ the modifications made in the original cascade model so that 
it applies to heavy nuclei. Some of the more fundamental con- 
cepts,' as they apply to the calculations, should be noted. Kirst 
of all, a nucleus having charge number Z and atomic mass number 
A consists of A nucleons originally traveling as a single particle. 
In the calculations it was assumed that each nucleon carried l/A 
of the total energy of the nucleus. 

A nucleus incident on an absorber was allowed to undergo 
its first interaction at a depth in the absorber selected random- 
ly from an exponential distribution which had a mean value equal 
to the mean free path of that nucleus in the absorber. The mean 
free path for a particular projectile nucleus y incident on a 
target nucleus x was calculated from the familiar expression for 

2 1 /"3 1 ^3 2 

the Interaction cross-section, = 7Tr^(Ay -^+ A^' - b) , 

where b is the so-called overlap parameter. For the calculations 
presented here, the value b = - 0.6 was used. Small changes ob- 
served in the overall cascade development by varying b are pre- 
sently being studied and will be discussed elsewhere. 

The various modes for generating nuclear fragments during 
the interaction of a heavy projectile nucleus with a heavy target 
nucleus are not yet well known. For the present calculations, the 
partial interaction cross-sections for heavy nuclei incident on 
hydrogen as calculated by Shapiro, et al . (1970) of the U.S. 

Naval Research Laboratory have been used to determine the proba- 
bilities for observing fragments with Z > 3. The probabilities 
for having no fragments with Z > 3 were extracted from the work 
of Cleghorn, et al. (1968). The number of alpha particles evolving 
from interactions of heavy nuclei were selected randomly from a 
Poisson distribution with mean value depending on the Z and A of 
the projectile nucleus as determined from emulsion measurements 
(Powell, et al. 1959). The numbers of protons and neutrons be- 
longing to the projectile (i.e. other than recoil nucleons) which 
leave the Interaction were determined from the conservation of 
nucleons after having first selected the type of fragment and 
the number' of alpha particles which left a particular interaction. 

It has been assumed that alpha particles and fragments 
with Z > 3 always left an interaction with an energy /nucleon 
equal to that of the primary nucleus at the instant of the inter- 
action. However, the nucleons were considered as having had 
the opportunity either (i) to participate in the interaction by 
exciting the target nucleus and producing secondary particles 
(plons), thereby losing part of their energy in the Interaction 
or (11) to leave the interaction with energy equal to the energy/ 
nucleon of the primary nucleus, i.e. not to participate in the 
interaction. 
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Each secondary heavy fragment was treated in the same 
manner as the Incident heavy nucleus (gen era tin;' additional 
heavy rragmi'nts, alpha particles, ar\d tuicleonr. 1. Tlu' alpha 
pai’ticies were allowed to initiate nucl ear-el ect.i’oniaiuud. 1 c 
cascades in the inannei’ described by Jones (KV/il, The nueleonr. 
were pei’mitted to undergo successive nuclear interactions, 
thereby producing independent nuclear-electromagnetic cascades. 

In other words, the cascade development initiated by an inci- 
dent heavy nucleus amounts to the superposition of the inde- 
pendent cascades being initiated at various depths in the ab- 
sorber by the nucleons comprising the original incident nucleus. 

3 . Results and Discussion. The average number of particles^ 

in the cascades initiated by and ^^i-'e nuclei incident on the 

top of the iron absorber (vertically, at the center) are shown 
in "Ig. 1. The results are given for the number of particles 
per nucleon of the incident nucleus for 30, 100 and joo GeV/nu- 
cleon total energy. For comparison, the cascade curves for 
primary protons are also given. As one might expect, the number 
of cascade particles per nucleon of the primary nucleus seems 
to be only weakly dependent on the atomic mass number of the 
nucleus. It appears that the cascade curve builds up and decays 
somewhat more slowly for Fe than for C, but the differences are 
certainly less pronounced than differences in the cascades for 
C and protons. The explanation for this feature is that the 
interaction length per nucleon is somewhat longer for heavier 
nuclei . 


Calculations are presently being carried out to study 
the effects of changes in the fragmentation parameters on the 
overall cascade development. The two extreme cases are, of 
course, (i) when a heavy nucleus loses only one nucleon in each 
consecutive interaction and (ii) when a nucleus is completely 
fragmented into its constituent nucleons in the first interaction. 
The most probable cascade development will be somewhere between 
these two extremes and should not be too different from the cal- 
culated results presented here. 

The vertical error bars shown in Fig. 1 should provide 
the reader with an idea of the fluctuations to be expected in 
the cascades at various stages in their development. Space 
does not permit a detailed presentation of these fluctuations. 
However, the Information given in Table I for the means and 
standard deviations of the sums of cascade particles per nucleon 
of the primary nucleus can be used to estimate the accuracy with 
which the cascade development can be used to make primary energy 


''' Here the term "n particles" is used to represent n times the 
ionization deposited in a scintillator by a 4 GeV muon tra- 
versing- (vertically) the scintillator. The energy loss rate 
of a 4 GeV muon is 17.4 % greater than that of minimum-ionizing 
muons. In other words, to obtain the correct multiple of the 
minimum-ionizing muon signal, one must multiply the given 
particle numbers by the factor 1.174. 
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Fig. 1. Mean numbers of cascade particles per nucleon of the 
primary nucleus as a function of absorber depth for 10, 30, 
and 300 GeV/Nucl protons (H), carbon nuclei (C), and iron 
nuclei (Fe) incident on an iron absorber. The top abscissa 
indicates the absorber depth in units of the proton interaction 
length (^)whlle the bottom abscissa indicates this depth in 
radiation lengths (r.l,). 
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Means and Standard Deviations (in parenthis) for the sums 
of cascade particles per nucleon of primary nucleus. Given as 
functions of Absorber Depth in Interaction Lengths {?[) and 
Primary Energy E in GeV. Particle Sampling is each 2 r • 1. 

The statistics cSlculated are given in the second column. 


E 

NUCLEUS 

2 ^ 

4;^ 

6;) 


lOi^ 

0 

50 

750 H 
100 C 
100 Pe 

221(131 ) 
206(65) 

204(41) 

360(98) 

348(43) 

352(25) 

397(76) 

394(25) 

399(11) 

405(67) 

405(22) 

410(12) 

407(66) 
408(21 ) 
412(13) 

100 

250 H 
100 C 
35 Fe 

679(401) 

605(243) 

569(169) 

1224(332) 

1175(175) 

1162(108) 

1411(213) 

1400(85) 

1399(47) 

1467(151) 

1465(48) 

1467(22) 

1481(142) 

l48l(4o) 

1482(17) 

300 

150 H 
40 C 
30 Fe 

1710(1153)5630(929) 

1495(711) 3425(482) 

1379(609) 3322(451) 

4355(583) 
4276(279) 
4237(201 ) 

4582(349) 

4586(143) 

4563(60) 

4636(304) 

4684(84) 

4655(35) 



2N [Particle^ 

Elg. 2. Distributions for the sums of cascade particles per 
nucleon (abscissa) for 100 GeV/Nucl carbon nuclei incident on 
an iron absorber. The distributions are for 5, 6, and 10 pro- 
ton interaction lengths (^) total absorber depth with the cas- 
cades being sampled each 2 radiation lengths (r.l.). 


„easure™ent= in .n^ron 

r/peL^or'tn; sSf fo/fncldent 100 OaV 0 nuclei 

for three different absorber depths. 
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It should be pointed out that the calculations presented 
herein are for the actual cascade particle numbers, l.e., no 
attempt has been made to take into account the inefficiency of 
the scintillators in recording the ionization of heavy Z parti- 
cles. If one were to consider this inefficiency, the cascade 
curves for C and Fe would be separated more in the early develop- 
ment, where a substantial number of heavy fragments contributes 
to the average cascade development. 
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